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Abstract 

The  calibrated  carrier-phase  GPS  technique  promises  time  transfer  between  remote  sites  with 
5-minute  precision  approaching  20  ps.  Uncalibrated  time  transfers  between  many  sites  are  routinely 
made  available  by  the  Earth  Orientation  Department  of  the  U.S.  Naval  Observatory  (USNO).  Although 
important  technical  issues  involving  long-term  accuracy  and  calibration  are  still  being  resolved,  it 
is  possible  to  generate  sample  time  and  frequency  scales  so  as  to  anticipate  the  full  power  of  the 
technique.  Using  minor  extensions  of  the  USNO's  postprocessed  extended-Percival  algorithm  SuperP 
[1],  a  carrier-phase-based  frequency  scale  is  generated  and  compared  to  local  USNO  timescales. 


INTRODUCTION 

The  U.S.  Naval  Observatory  (USNO)  is  actively  pursuing  the  concept  of  a  “Distributed  Master 
Clock,”  which  would  use  GPS  carrier-phase  data  to  provide  on-line  timing  ensembling  and 
synchronization  of  many  remote  sites.  It  is  anticipated  that  easily  incorporated  improvements 
to  the  services  currently  provided  by  the  Analysis  Centers  of  the  International  GPS  Service 
(IGS)  would  prove  adequate  for  this  purpose.  Among  the  required  improvements  would  be  the 
use  of  calibrated  and  calibratable  GPS  receivers  [2]  and  the  elimination  of  solution  day-boundary 
discontinuities.  Although  these  improvements  are  all  necessary  or  important  for  time  transfer, 
they  do  not  obviate  frequency  transfer,  or  the  generation  of  frequency  scales.  We  of  course 
do  not  wish  to  trivialize  the  difficulties  or  the  practical  importance  of  being  able  to  integrate  a 
frequency  scale  to  a  timescale;  this  work  should  therefore  be  considered  a  progress  report  on 
efforts  to  assess  the  promise  of  the  carrier-phase  technique. 


DATA  REDUCTION 

The  input  data  for  this  analysis  were  obtained  from  the  daily  reductions  of  carrier-phase  data 
routinely  produced  by  the  USNO’s  Earth  Orientation  (EO)  Department,  and  initially  included 
47  sites  whose  timing  reference  was  listed  as  a  maser  in  their  site  logs,  as  listed  by  the  IGS 
with  URL  http://igscb.jpl.nasa.gov/network/site.  The  data  were  taken  from  MJD  51062  through 
51384  (6  Sep  98-25  Jul  99).  Timing  data  from  the  IGS  sites  designated  USNO  (located  at  our 
Washington,  DC,  site)  and  AMC1  and  AMC2  (located  at  the  USNO’s  Alternate  Master  Clock 
at  Schriever  AFB  in  Colorado)  were  automatically  corrected  for  the  TurboRogue  receivers’ 
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internal  delay  jumps  using  the  procedures  outlined  in  [3],  but  this  made  little  difference  as  the 
timing  data  were  then  replaced  by  their  first  differences  (frequencies)  that  were  then  edited  for 
outliers.  Day-boundary  discontinuities  were  avoided  by  excluding  frequencies  generated  from 
data  that  spanned  those  boundaries,  and  only  reductions  in  which  the  USNO  receiver  was  used 
as  a  reference  were  utilized.  Outlier-removal  algorithms  eliminated  all  single-point  fluctuations 
in  excess  of  15  “deviations,”  where  one  “deviation”  is  defined  to  be  half  the  difference  between 
the  frequency  exceeded  by  84%  of  the  data  minus  the  frequency  exceeded  by  16%  of  the  data 
(so  that  the  “deviation”  is  to  the  standard  deviation  as  the  median  is  to  the  average). 

The  clock  data,  in  the  form  of  estimated  and  edited  frequency  differences  between  the  USNO 
and  site  receivers,  were  then  input  into  the  USNO  program  SuperP  for  clock  characterization  [1]. 
This  program  applies  a  postprocessed  Percival-type  algorithm  which  computes  the  instrumental 
frequency  bias  and  drift  of  all  clocks  relative  to  one  fiducial  clock  via  least  squares.  Data  from 
clocks  whose  frequency  or  rate  seemed  to  shift  suddenly  to  a  new  value,  relative  to  the  USNO 
reference,  were  treated  as  if  from  a  different  clock.  The  time  reference  for  the  USNO  receiver 
is  Master  Clock  #3  (MC3),  and  is  referred  to  as  the  “USNO  reference”  because  it  was  only  used 
as  a  reference,  and  has  no  influence  on  the  clock  characterization  or  the  subsequent  averaging 
for  the  generation  of  a  frequency  scale.  The  absolute  frequency  and  frequency  drift  biases  of 
the  fiducial  clock  cannot  be  determined  from  clock  differences;  however,  for  the  purpose  of 
data  editing  they  were  chosen  so  as  to  keep  the  average  of  individual  clock  frequencies  and 
drifts  close  to  the  USNO  reference.  Because  these  assumptions  and  this  procedure  result  in  an 
undetermined  overall  rate  and  drift  shift  in  the  frequency  scales  produced,  the  standard  “<7y” 
analysis  has  been  supplemented  with  the  “az”  statistic  [4,5],  which  is  independent  of  both  rate 
and  drift.  Both  these  measures  were  generalized  slightly,  as  appropriate,  to  allow  for  the  fact 
that  the  data  are  unevenly  spaced  and  expressed  in  terms  of  frequency  instead  of  phase. 

It  was  found  that  outlier  removal  on  the  basis  of  individual  5-minute  frequency  points  was 
insufficient  to  identify  all  incorrect  data,  as  there  were  times  when  sequences  of  obviously 
erroneous  frequencies  would  be  revealed  only  after  the  data  were  averaged  into  hourly  or  half¬ 
day  bins.  In  some  cases,  the  site  logs  showed  that  these  variations  were  related  to  maser  tuning 
or  repair;  however,  at  other  times  there  was  no  such  indication.  Although  we  anticipate  that 
at  some  point  the  IGS  will  stipulate  that  clock  records  be  systematically  maintained  and  made 
publicly  available,  there  will  always  be  unexpected  frequency  excursions  related  to  hardware  or 
software  failures.  For  this  work,  outlier  removal  was  completed  by  a  manual  iterative  process 
in  which  SuperP-corrected  clock  data  were  inspected,  and  then  new  clock  characterizations 
generated  on  edited  data.  Clocks  whose  frequencies  showed  significantly  higher  variance  or 
extreme  non-white  behavior  were  excluded  completely,  so  that  in  the  end  only  data  that  were 
used  came  from  the  following  14  IGS  sites:  ALGO,  AMC,  DRAO,  FAIR,  GODE,  HOB2, 
HRAO,  KOKB,  MATE,  MKEA,  NLIB,  NRC1,  PIE1,  and  USNO.  Figure  1  shows  how  many 
sites  contributed  to  the  final  frequency  scales  at  any  one  time. 

One  potentially  serious  issue  is  that  for  some  sites  we  were  not  able  to  obtain  information 
concerning  clock  steering.  Although  the  SuperP  algorithm’s  implementation  should  remove 
the  effects  of  large,  sudden  changes,  any  clocks  being  smoothly  steered  to  UTC  would  simply 
appear  to  be  very  good  clocks  over  the  long  term.  However,  since  the  USNO  unsteered 
maser  mean  timescale  (UNSTMM)  is  probably  more  frequency  stable  than  UTC  over  the  time 
range  of  interest  (due  to  the  noise  of  common-view  time  transfer  being  folded  into  UTC),  any 
unmodeled  steering  of  IGS  site  clocks  would  tend  to  weaken  the  conclusions  of  this  paper,  and 
proper  allowance  for  such  steering  would  strengthen  them. 

For  clock  characterization,  fits  were  made  to  frequency  differences  between  each  pair  of  clocks, 
and  points  were  weighted  by  the  inverse  of  their  squared  formal  errors,  as  computed  statistically 
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assuming  the  error  values  output  by  GIPSY  were  uncorrelated.  It  was  found  that  essentially 
the  same  results  were  achieved  using  unit  weights  for  clock  characterization. 


THE  FREQUENCY  SCALE 

Once  the  iterative  process  of  clock  characterization  was  complete,  detrended  clock  frequencies 
were  averaged  into  3-hour  bins.  This  was  necessary  to  make  the  subsequent  average  less 
sensitive  to  the  intermittent  loss  of  data  that  have  a  bias.  Not  surprisingly,  for  data  massaged 
so  as  to  maximize  their  “whiteness,”  it  was  found  that  weighting  data  by  their  inverse  variances 
produced  a  closer  match  between  the  individual  clock  frequencies  and  the  scale  than  did 
unit-weighting.  However,  it  also  produced  a  closer  match  between  the  frequency  scale  and  the 
UNSTMM,  which  played  no  role  in  the  frequency  scale  generation  (nor  did  any  of  the  USNO 
masers).  Figure  2  shows  that  the  frequency  scale  was  closer  to  our  UNSTMM  than  to  USNO 
MC3  (steered  to  the  USNO  Master  Clock,  which  in  turn  is  steered  to  UTC),  and  from  this  we 
conclude  that  the  frequency  scale  derived  from  IGS  data  is  meaningful. 

Figure  3  displays  the  ay  plots  associated  with  the  scales  in  Figure  2.  A  frequency  drift  was 
removed  from  these  plots,  which  forces  each  curve  to  negative  infinity  for  large  r. 

Figure  4  displays  az  plots  with  the  scales  of  Figure  2.  az  can  be  considered  a  generalization  of 
the  Hadamard  variance  for  unequally  spaced  data. 

Figure  5  compares  ay  plots  of  individual  IGS  sites,  using  the  IGSFS  as  a  reference,  with  the  ay 
plots  of  the  local  USNO  masers,  using  the  UNSTMM  as  a  reference,  and  Figure  6  uses  the  az 
statistic  for  the  same  data. 

FREQUENCY  TRANSFER  TO  THE  USNO  ALTERNATE  MAS¬ 
TER  CLOCK  (AMC) 

Data  from  the  USNO  AMC  were  not  included  in  the  frequency  scale  analysis  because  its 
master  clock  is  steered  to  UTC(USNO)  and  we  wished  to  avoid  even  the  suspicion  of  adding 
pre-correlated  data.  However,  the  AMC  is  uniquely  qualified  to  serve  as  a  time-transfer  test 
bed  due  to  its  large  ensemble  of  clocks  and  the  hourly  time  transfers  with  USNO-DC  via 
calibrated  Two-Way  Satellite  Time  Transfer  (TWSTT). 

Figures  7  and  8  display  the  ay  and  az  plots  of  the  differences  between  the  UNSTMM  (unsteered 
average  of  up  to  10  masers  at  USNO-DC)  with  the  unsteered  average  of  the  two  masers  at 
the  USNO  AMC  when  frequency  (or  time)  transfer  is  achieved  via  carrier-phase,  P-code 
common-view  GPS,  and  TWSTT.  Common-view  data  are  from  48-hour  linear  fits  to  13-minute 
single-channel,  dual-frequency  observations.  The  downconverter  in  the  AMC  GPS  receiver’s 
antenna  was  improved  subsequent  to  this,  and  more  recent  data  indicate  a  performance  better 
by  about  .25  in  logarithmic  units. 


OCCASIONALLY  IGNORED  TIME-TRANSFER  ISSUES 

The  primary  problem  in  implementing  carrier-phase  time  transfer  is  for  calibrated  and  calibrat- 
able  receive  systems,  which  are  both  environmentally  protected  and  environmentally  insensitive, 
such  as  the  Geodetic  Time  Transfer  receivers  (GeTT,  designated  USNB  in  the  data  reduction) 
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[6].  Herein  we  discuss  two  other  considerations:  carrier-phase  data-analysis  noise,  and  errors 
in  ensembling  local  clocks. 

“Analysis  noise”  is  a  term  loosely  used  to  describe  the  fact  that  different  researchers  or  op¬ 
erational  centers  can  compute  different  time-transfer  results  between  two  sites  despite  using 
identical  data  from  those  sites.  This  can  be  due  to  differences  in  software,  modelling,  pa¬ 
rameterizing,  editing,  or  use  of  data  from  other  sites  in  a  combined  solution.  In  order  to 
provide  a  crude  estimate  of  analysis  noise,  data  reduced  by  the  USNO  EO  Department  using 
the  GIPSY  software  package  were  compared  to  data  reduced  by  the  CODE  analysis  center  in 
Bern,  Switzerland,  using  the  Bernese  software  package.  Figure  9  compares  double-differenced 
data  from  the  short-baseline  USNO  and  USNB  receivers  with  data  from  the  2000-km  baseline 
AMC  vs.  USNO  and  AMC  vs.  USNB.  Large  excursions  are  typically  associated  with  high 
formal  errors  from  the  solutions  and  have  been  excluded  in  the  figure.  Excluding  those  outliers, 
we  estimate  about  1  ns  for  the  peak-to-peak  analysis  noise  in  GPS  carrier  phase. 

The  largest  error  in  the  USNO-DC  measurements  is  due  to  the  mostly  underground  time 
transfer  between  our  two  buildings,  which  is  realized  via  the  phase  of  the  clock  5  MHz  signals 
(Table  1).  Operational  data  are  transferred  to  Building  78  from  Building  52  via  fiber  optics 
using  an  LED  light  source.  Two  additional  signals  are  also  sent  via  parallel  means:  MC3’s 
signals  (in  Bldg.  52)  are  sent  through  a  different  set  of  fiber  optics  using  lasers,  and  signals 
from  one  maser  (NAV2)  are  sent  to  Bldg.  78  via  a  coaxial  cable.  In  the  reverse  direction, 
lasers  are  used  to  send  the  signal  from  the  USNO  MC  (MC2,  in  Bldg.  78)  to  Bldg.  52.  The 
laser  signals  are  also  split  at  each  end  so  as  to  be  measured  via  multiple  measurement  systems, 
which  agree  to  about  100  ps,  peak  to  peak.  Since  we  also  have  TWSTT  and  carrier-phase 
connections  between  the  two  buildings,  the  USNO-DC  is  an  excellent  test  bed  for  short-baseline 
time  transfer. 

By  differencing  the  values  of  MC2-MC3  measured  by  each  signal  going  to  Bldg.  78  to  create 
double-differences,  it  is  possible  to  measure  the  relative  change  in  their  delay.  It  is  also  possible 
to  measure  the  variations  of  the  total  round-trip  delay,  by  differencing  MC2-MC3  from  the 
laser-based  signals  going  to  Bldg.  78  from  those  going  from  Bldg.  78.  Figure  10  shows  these 
differences  along  with  the  measured  temperatures  along  the  path,  and  Figure  1 1  compares  the 
same  double-differenced  delays  with  the  measured  relative  humidities. 

The  round-trip  differences  provide  the  best  measure  of  the  change  in  delay  with  temperature, 
which  would  tend  to  cancel  for  signals  going  in  the  same  direction.  The  outdoor  temperature 
appears  to  be  the  only  significant  contributor;  some  apparent  correlations  with  relative  humidity 
are  probably  due  to  associated  temperature  variations.  The  temperature  of  the  room  in  which 
the  time  distribution  electronics  and  measurement  systems  are  kept  appears  to  be  significant 
only  at  times  of  >10°F  environmental  variations  due  to  control  equipment  failure  (which  are 
not  shown  in  the  figure).  A  variable  diurnal  cycle  of  roughly  70-150  ps  peak-to-peak  amplitude 
is  correlated  with  the  time  of  the  exterior  daily  temperature  variations,  but  not  the  amplitude 
(Figure  12).  Comparisons  of  the  transmitted  signals  with  those  of  unsteered  clocks  at  each  site 
show  that  roughly  two-thirds  of  the  round-trip  delay  is  in  the  line  from  Bldg.  78  to  Bldg.  52, 
and  one-third  in  the  lines  from  Bldg.  52  to  Bldg.  78. 

The  observed  variation  is  too  large  to  be  explained  by  thermal  expansion  of  the  glass:  assuming 
a  cable  length  of  160  m  underground  and  70  m  within  buildings,  a  glass  index  of  refraction  of 
1.4,  a  refraction  angle  correction  of  1.2,  and  a  thermal  expansion  coefficient  of  8.5  x  10~6/°C,  one 
expects  a  round-trip  delay  variation  of  15  ps/°C  underground  and  7  ps/°C  within  the  buildings. 
The  observed  underground  seasonal  variations  of  22° C.  are  an  order  of  magnitude  too  small  to 
explain  the  observed  3  ns  seasonal  variations  in  round-trip  delay.  Although  these  calculations 
do  not  include  the  dielectric  constant,  they  appear  to  be  consistent  with  previously  reported 
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values  of  6-7  ppm/°C  [7-9].  In  the  context  of  this  paper  our  inability  to  explain  the  variation  is 
not  important,  because  we  have  recently  received  and  plan  to  install  temperature-compensated 
coaxial  cable  and  fiber  optics,  and  also  because  the  observed  associated  single-direction  delay 
variations  are  less  than  1  ns  peak  to  peak,  and  therefore  adequate  for  the  purpose  of  evaluating 
carrier  phase  to  that  level.  Figure  13  shows  that  carrier-phase  timing  differences  measured 
using  receivers  USNO  (in  Bldg.  52)  and  USNB  (the  GeTT,  in  Bldg.  78)  show  much  larger 
variations,  of  about  3  ns  peak  to  peak.  This  is  probably  all  due  to  errors  in  the  carrier  phase 
technique  as  currently  practiced,  and  it  is  interesting  that  these  errors  are  similar  to  those  in 
the  TWSTT  link  between  the  two  buildings. 


CONCLUSIONS 

We  conclude  that  it  is  currently  possible  to  generate  an  excellent  frequency  scale  using  GPS 
carrier-phase  data,  and  anticipate  no  significant  problems  in  the  generation  of  a  timescale, 
once  known  technology  is  actually  installed  at  a  sufficient  number  of  sites.  Analysis  noise 
is  subnanosecond  rms,  and  hardware  limitations  may  be  1  ns  rms,  ignoring  currently  serious 
hardware  calibration  considerations  in  our  equipment. 
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Table  1.  USNO-DC  Link  Between  Buildings 


Bldg.  52  < - >  Bldg.  78 

MC3  =  IGS  ref  MC2  =  UTC(USNO) 


I.  5  MHz  signals  sent  from  Bldg.  52  to  Bldg.  78 

A.  All  clocks  via  LED -based  fiber  optics 

B.  MC3  via  laser-based  fiber  optics  (B  -  A  <  1  ns) 

C.  One  maser  CNAV2')  via  RG223  coaxial  cable  (C  -  A  <  2  nsl 

D.  TWSTT  (D'-  A  <3  ns) 

E.  Carrier-phase  (USNO  -  USNB)  (E  -  A  <  3  ns) 

F.  Portable  cesium  1-pps  (F  -  A  <  .5  ns) 

(Values  in  parentheses  estimate  peak-to-peak 
differences  between  each  technique  and  the  first). 

II.  5  MHz  signal  sent  from  Bldg.  78  to  Bldg.  52 
via  laser-based  fiber  optics 

Round  Trip-Delay  Variations: 

(MC2  -  MC3)bm3.52  -  (MC2  -  MC3)BWg.78 

Round-Trip  Delay  <  3  ns  (seasonal  p/p) 

<  150  ps  (diurnal  p/p) 

C.  Temperature-compensated  coaxial  cables  (to  be  installed) 
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Figure  1.  Histogram  of  number  of  sites  contributing  to  frequency  scale  as  function  of  time. 


Figure  2.  Frequency  scales,  in  units  of  parts  in  1015  (fs/s).  The  top  plot  is  the  unsteered  USNO 
maser  mean  minus  the  carrier-phase  based  scale;  middle  is  the  unsteered  maser  mean  minus 
MC3;  and  bottom  is  the  carrier-phase  frequency  scale  minus  MC3. 
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oy(t)  after  frequency  drift  removed 


Figure  3.  ay  plots  of  parabola-removed  frequency  scales  in  Figure  2.  Symbol  1  is  the  carrier- 
phase  frequency  scale  minus  MC3;  symbol  2  denotes  the  unsteered  maser  mean  minus  Master 
Clock  3;  and  symbol  3  denotes  the  unsteered  maser  mean  minus  carrier-phase  frequency  scale. 


Figure  4.  ax  plots  of  frequency  scale  differences;  symbols  have  same  significance  as  in  previous 
figure. 
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oy(t)  of  frequency  drift  removed  IGS  clocks  -  tGS  ensemble 


tog(t  seconds) 

Figure  5.  Top:  ay  of  each  detrended  IGS  maser,  using  the  IGS  frequency  scale  as  a  reference. 
Middle:  ay  of  each  edited  and  detrended  USNO  maser,  using  the  USNO  unsteered  maser- 
mean  (UNSTMM)  as  a  reference.  Bottom:  uy  for  uncorrected  USNO  maser  data,  using  the 
UNSTMM  as  a  reference. 


Iog(x  seconds) 


Figure  7.  av  plots  of  parabola-removed  timing  differences  between  the  USNO  unsteered  maser 
mean  and  the  average  of  the  two  AMC  masers.  Symbols  1  and  2  use  carrier-phase  via  IGS 
receivers  USNO  or  USNB;  symbol  3  uses  dual-frequency  P-code  GPS  common  view;  and  symbol 
4  depicts  TWSTT. 


oz(r)  AMC  -  USNO  unsteered  maser  means 


Figure  8.  As  in  previous  figure,  except  <rz  is  used  instead  of  cry. 


Analysis  Noise:  USNO-EO  vs  CODE  Bernese 


Fiber  delays  vs  temperature 


Figure  9.  Analysis  noise, 
as  measured  by  differ¬ 
ences  in  solutions  gener¬ 
ated  by  USNO  and 
CODE.  Top  is  for  base¬ 
line  USNO-AMC1,  and 
bottom  is  for  USNB- 
USNO.  Very  large  excur¬ 
sions,  not  shown,  are 
associated  with  high 
formal  errors  computed 
by  the  solution. 


Figure  10.  Timing  differ¬ 
ences  between  Bldgs.  78 
and  52  in  USNO-DC  in 
the  form  of  MC2-MC3 
measurements  subtracted 
from  MC2-MC3  meas¬ 
urements  using  an 
LED -based  fiber-optic 
system.  The  upper  three 
plots,  from  top  to 
bottom:  uncalibrated 
TWSTT,  underground  co¬ 
axial  cable,  and  laser- 
based  fiber  optics. 

The  coaxial-LED  double¬ 
difference  data  in  the 
second  plot  extend  back 
through  4  years  to  MJD 
49800,  and  show  a  fairly 
constant  peak-to-peak 
seasonal  variation  of  2 
ns.  The  fourth  plot  is 
the  round-trip  delay 
measured  by  two  laser- 
links.  The  lower  four 
plots,  from  top  to  bot¬ 
tom:  temperatures  in  the 
Bldg.  78  measurement 
room,  near  an  under¬ 
ground  tunnel  junction, 
outdoors,  and  in  the 
Bldg.  52  measurement 
room. 
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Fiber  delays  vs  relative  humidity 


MJD 


Figure  11.  The  upper  four  plots  are  the  same  timing  differences  as  in  Figure  10.  The  lower 
three  plots,  from  top  to  bottom:  relative  humidity  of  the  Bldg.  78  measurement  room,  outdoors, 
and  the  Bldg.  52  measurement  room. 
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Figure  12.  Round-trip  path  variations  between  USNO  Buildings  (top)  and  the  exterior  tern- 
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nanoseconds 


USNO-USNB  common  clock  corrected  by  CODE 
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Formal  Errors  of  USNO  EOG  solution 
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Figure  13.  Upper  plot  is  double  difference  of  timing  differences  between  USNO  Bldgs.  52 
and  78  measured  by  carrier  phase  and  reduced  by  CODE/Bernese  vs.  those  measured  using 
LED-based  fiber  optics.  Constants  were  removed  to  handle  GPS  receiver  jumps  just  before 
MJDs  51038,  51051,  and  51226.  Lower  plot  is  the  formal  error  of  the  GIPSY  computations. 
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Questions  and  Answers 


PETER  WOLF  (BIPM):  This  is  just  a  detail.  You  said  that  at  at  some  point  that  you  cheated 
for  phases,  you  called  it;  you  smoothed  over  date  boundaries.  That  means  you  just  adjusted 


days 


W/k^t/*n1lT7 

iiiw  j  ump,  i/aaiva.ii y . 


DEMETRIOS  MATSAKIS  (USNO):  What  I  actually  did  was  convert  everything  to  frequencies 
and  then  throw  out  the  frequency  that  crossed  the  day  boundary. 


GERARD  PETIT  (BIPM):  By  doing  so,  don’t  you  extend  the  short-term  stability  of  the  maser 
to  longer  term?  By  removing  the  differences  at  the  boundary? 

MATSAKIS:  Well,  that  is  just  a  difference  between  two  5-minute  points.  So  basically.  I’m 
assuming  that  the  maser  is  constant  over  those  10  minutes,  constant  in  rate  and  drift.  If  a 
maser  did  something  very  funny  right  there  and  it  just  happened  to  correspond  to  midnight, 
that  would  certainly  be  hidden. 


But  actually  the  Percival  algorithm  washes  all  those  problems  away  anyway. 

DAVID  ALLAN  (Allan’s  Time):  Had  you  looked  at  JPL -calibrated  round-trip  fiber  optic? 
They’re  doing,  I  think,  30  picoseconds  at  the  antenna  with  respect  to  their  masers,  using  that 
technique  for  comparing  the  two. 

MATSAKIS:  I’m  aware  of  them.  I’m  not  aware  of  the  bill.  It’s  how  much  do  we  want  to  spend. 
And  as  we  get  more  into  carrier  phase,  the  kind  of  money  those  kinds  of  systems  involve  gets 
less  and  less,  figuratively  speaking. 

We  also,  of  course,  have  plans  in  5  or  10  years  to  have  a  new  building  which  will  house  all  our 
clocks.  And  all  the  fiber-optic  questions  will  go  away.  Although  it’s  an  incredible  annoyance 
for  us  to  deal  with  fiber  optics,  the  reason  I  brought  it  up  was  because  I  was  more  interested  at 
looking  at  what  short-baseline  carrier  phase  does.  And  to  show  that  the  fiber  optics,  although 
they  were  a  problem,  were  less  of  a  problem  than  the  carrier  phase  was  itself. 
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